Cluster Probability in Bootstrap Percolation 
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We develop a recursive formula for the probability of a fc-cluster in bootstrap percolation. 
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I. INTRODUCTION 

The model of bootstrap percolationi (or fc-core perco- 
lation, as we will refer to it) has a connection to jam- 
ming transitions, including the glass transition^ Many 
features of dynamical arrest in glassforming liquids are 
captured by kinetically-constrained spin models,'^ which 
map onto fc-core percolation and its variants. In addi- 
tion, Schwarz, et al^ have argued that there is an anal- 
ogy (but not a strict mapping) between /c-core percola- 
tion and the jamming transition of frictionless granular 
packings^'^, similar to the analogy between fc-core per- 
colation and rigidity percolation.^ In the mean-field, or 
infinite-dimensional limit, the fc-core percolation transi- 
tion was shown to be peculiar in Ref . ^ and subsequently 
in Ref. 0. The transition is discontinuousfi*^ but also 
exhibits power-law scalingiiS and a diverging susceptibil- 
ity and length scaled. This unusual type of transition 
has been called a random first-order phase transition^ 
What makes the mean-field fc-core percolation transition 
particularly interesting is that its exponents are identical 
to the mean-field exponents of several models that have 
been proposed for jamming transitions, ^'^•^^■^^'-'^^ raising 
the possibility that fc-core percolation may be in the same 
universality class as these models. 

For jamming transitions, the random first-order nature 
of the transition does not appear to be an artifact of the 
mean-field approximation. Numerical simulations^'^ and 
scaling argumentsi^*i^ suggest that the jamming transi- 
tion of frictionless granular packings (Point J) is discon- 
tinuous with a diverging length scale in two and three 
dimensions. It has also been argued that the finite- 
dimensional glass transition, if it exists as a true phase 
transition, has a mixed character, with a discontinuity in 
the infinite-time limit of the dynamical structure factor 
and a diverging dynamical length scaleiiSiii This raises 
the question of the nature of the fc-core percolation tran- 
sition in finite dimensions d. Expansions in l/d indicatai^ 
that that the mixed nature of the transition should per- 
sist for a range of large but finite d. In two dimensions, 
specific variants of fc-core percolation have been identified 
that are rigorously known to exhibit discontinuous tran- 
sitions with diverging length scales, both at p = 1 (full 
occupancy and p < 1 (partial occupancy) li^ Numer- 
ical studies of other variants^ also point to the possibility 
of mixed transitions in finite dimensions at p < 1. How- 
ever, it is also known that some realizations of the model, 
such as fc = 3 on the two-dimensional triangular lattice. 



exhibit continuous transitions Thus, in finite dimen- 
sions, it appears that different versions of fc-core models 
can exhibit at least three different types of transitions: 
continuous transitions, mixed transitions at p = 1 and 
mixed transitions at p < 1. 

In order to sort through this confusing state of affairs, 
it would be useful to map fc-core percolation onto a spin 
model, which could then be studied by a wide variety of 
methods. Such an approach proved extremely powerful 
for ordinary percolation. A necessary first step towards 
this goal is to be able to count fc-core clusters with the 
correct probabilities. In this paper, we show how this 
can be done by deriving a relation for the probability 
distribution of clusters. 

The fc-core ensemble is defined as follows. We start 
with the ensemble of 2^ percolation configurations on a 
finite lattice of TV bonds. Each percolation configuration 
Cp has probability 



P 



(1) 



where riocc(Cp) and ?T.vac(Cp) are, respectively, the number 
of occupied and vacant bonds in the percolation config- 
uration Cp. A fc-core configuration Ck is obtained from 
Cp by removing all bonds which intersect a site to which 
there are less than fc remaining bonds. This process is 
recursively continued because as bonds are so removed, 
they may cause more sites to have less than fc neighbors. 
After this culling, the remaining fc-cluster is one in which 
all sites have at least fc occupied bonds. The weight of 
the configuration Ck is the sum of the percolation proba- 
bilities over all the percolation configurations which gave 
rise to Ck- 

In the percolation problem one can easily establish that 
the probability of finding an isolated cluster— Tp is given 

by 



i(rp) 



t(rp 



(2) 



where n{Tp) is the number of occupied bonds in the clus- 
ter Tp and t{Tp) is the number of perimeter bonds which 
must be unoccupied in order to define the limit of the 
percolation cluster Tp. One would like to write a similar 
formula for Pk{T), the probability of a fc-cluster Ffc, but 
it is obvious that this is not quite so trivial because of the 
complications of the culling process. It is the purpose of 
this paper to develop such a formula. 
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II. FORMULATION 

To proceed we introduce the function er(Cp) which is 
defined as unity if the percolation configuration Cp, when 
cuUed, gives rise to the rooted fc-cluster F. A rooted 
cluster includes the site at the origin and two rooted 
clusters are distinct unless their list of occupied bonds 
coincides exactly. The final culled configuration must 
contain the rooted fc-cluster F, but can also contain an 
arbitrary number of /c-clusters disconnected from F. If 
culling of Cp does not produce a configuration contain- 
ing the rooted cluster F, then er(Cp) is zero. With this 
definition is seems obvious that 



Note that the sum over fc-clusters F' does not include 
F' = F. Equation iQ is the principal result of this pa- 
per. On a qualitative level we can see the following: 
for percolation Eq. l^Jl gives what could be called a 
perimeter renormalization (1 — p)"'*^^\ whereas the sub- 
traction terms in Eq. iQ) give a much weaker "culling 
renormalization." In the extreme case of lattice animals, 
there is no perimeter renormalization at all. In that case 



P{C) 



P 



(where p is interpreted to be the bond 



Ffc(F) - ^er(Cp)P,(Cp) . 



(3) 



This relation expresses the fact that the cluster F inher- 
its the total percolation probability of all configurations 
which contain the cluster F after culling. Note that a 
fixed F inherits probability from not only single clusters 
which, when culled, yield F, but also from the vast array 
of configurations which, when culled, yield F together 
with arbitrary disconnected fc-clusters. Now it is clear 
that the sum over Cp can be restricted to configurations 
which either are equal to F or include it as a proper sub- 
set. So we write 

Pfc(F) = J2 er{Cp)Pp{Cp) 

= E (l-[l-er{Cp)])Pp{Cp) . (4) 



fugacity, usually denoted K.). The continuous aspect of 
fc-core percolation haa^ the same anomalous correlation 
length exponent v — 1/4 as does lattice animals in infi- 
nite dimensions p22i A heuristic explanation of this is that 
the culling renormalization is probably closer to the ab- 
sence of a perimeter renormalization in lattice animals 
than to the perimeter renormalization of percolation. 

Now we briefiy explore some consequences of the above 
result. We iterate Eq. Q to get 
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But 



Cp '.Cp 



E 

Crj - Cqr) 



where, for simplicity, we do not write out (in the last 
line) the inclusions. Because we are dealing with a finite 
lattice, the index m is bounded by a maximum value, 
denoted mmax- 

Next we consider the calculation of the connectedness 
susceptibility, %, which can be written as 



Pp {Cp 



(5) 



X 



(9) 



because after we occupy the bonds of F (with probability Thus we write 
^"(1")) we sum over all states of all the bonds not in F 
(which are occupied with probability p and vacant with 
probability 1 — p). Thus 



p 



«(r) 



E 

Cp 'Cp ^ r 



[1 - er{Cp)]Pp{Cp) . (6) 



E(-i)" E 

{rfc(m)} 



m— 1 



(10) 



The result of Eq. |(5J for the probability of a percola- 
tion cluster can be calculated in this same style. Note 
the meaning of the factor [1 — er(Cp)]: it is an indicator 
function such that culling does not lead to the fc-cluster 
F. But since Cp contains F but does not cull to F, it must 
be that the configuration Cp leads to a cluster F> which 
includes F as a proper subset. So we have 



By relabeling clusters, this can be rearranged to combine 
all terms with the same factor ^"(i'), in which case one 
has 



E^^ 



(r^) 



(1)m2 



rS,^'<r 



P,(F) = ^ J2 er'{Cp)Pp{Cp) 

r':r'>rCp:Cp>r' 

r,«(r) _ 
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(11) 



3 



This result no longer requires that Af be finite because 
for any finite cluster F^;, the number of terms coming 
from included subclusters in the sum inside the square 
brackets is independent of system size. Thus this for- 
mula presumably can be used as long as p < pk and it 
may yield a possible route to an analytic approach to the 
calculation of x for finite dimensional lattices. 

The result Eq. can be generalized to any observ- 

able G(rfc) that can be defined for a k-cluster F^. In that 
case, the ensemble-averaged value of G is 

where the cumulant function Gc is 



... + (-ir E gcfI")) 

r<'"'<r(— i'...<r<i)<r, 
+ ... (13) 

or more simply, 

G,(F,) = G(Ffe) - Gcir'i^). (14) 

This cumulant subtraction for the fc-core percolation 
problem is much more difficult to handle than it is for 
ordinary percolation. 

We thank Vincenzo Vitelli for stimulating discussions 
and a critical reading of the manuscript. This work was 
supported by NSF-DMR-0605044 (AJL). 
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